Unconventional continuous phase transition in a three dimensional dimer model 
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Phase transitions occupy a central role in physics, due both to their experimental ubiquity and 
their fundamental conceptual importance. The explanation of universality at phase transitions 
was the great success of the theory formulated by Ginzburg and Landau, and extended through 
the renormalization group by Wilson. However, recent theoretical suggestions have challenged this 
point of view in certain situations. In this Letter we report the first large-scale simulations of a 
three-dimensional model proposed to be a candidate for requiring a description beyond the Landau- 
Ginzburg- Wilson framework: we study the phase transition from the dimer crystal to the Coulomb 
phase in the cubic dimer model. Our numerical results strongly indicate that the transition is con- 
tinuous and are compatible with a tricritical universality class, at variance with previous proposals. 
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The Landau-Ginzburg- Wilson (LGW) theory of phase 
transitions Q has been a remarkably powerful approach 
to study critical phenomena, both in classical and quan- 
tum systems. The basic assumption is that each phase 
is characterized by its symmetries and a local order pa- 
rameter for the broken symmetries (in case of an ordered 
phase). To describe a phase transition, this approach 
amounts to expanding the free energy of the system in 
powers of the order parameter(s) describing the ordered 
phase(s). Exceptions are however known in 2d classical 
or equivalently Id quantum systems where stable critical 
phases (with power-law correlations but no broken sym- 
metry) exist, the low temperature phase of the 2d XY 
model being one of the most famous examples. The asso- 
ciated Kosterlitz-Thouless (KT) phase transition and 
the role of topological defects have become very fruitful 
concepts in the statistical physics of 2d phases not de- 
scribed by a simple order parameter in the LGW sense. 

The stimulating idea that "non-LGW" continuous 
phase transitions could also occur in higher dimensions 
has recently been proposed in the context of quantum 
magnetism There, the possibility of a generic {i.e. 
not requiring fine tuning) continuous transition between 
two phases with different symmetry breakings (Nccl and 
Valence Bond Solid states) was pointed out, in contrast to 
the LGW prediction of cither an intermediate phase with 
none of the two orders, or a coexistence region, or a direct 
first order transition. So far, simulations on all candidate 
quantum models rather see a first order phase tran- 
sition, a possibility which can never be discarded for a 
given microscopic model (see, however, recent claim of a 
continuous transition [^). 

From these perspectives, classical dimers at close pack- 
ing on simple hypercubic lattices are particularly inter- 
esting as they are too constrained to form a liquid with 
a finite correlation length, even at infinite temperature: 
instead of decaying exponentially, the dimer-dimer cor- 



relations are algebraic, both in 2d [6j and 3d |9|. In 
2d, the transition from such a critical state to a bro- 
ken symmetry phase (dimer crystal) has been studied in 
Ref. |y. In this Letter we study the analogous transi- 
tion in a three-dimensional classical dimer model. Both 
in 2d and 3d, a naive LGW expansion in terms of the 
low-tcmpcrature order parameter completely misses the 
critical nature of the high-temperature phase and thus 
cannot describe correctly the transition. Whereas in 2d 
the situation is well understood in terms of a KT transi- 
tion 1^ , our numerical simulations on the 3d model show 
the existence of a single continuous second-order phase 
transition between the algebraic liquid at high tempera- 
ture (so-called Coulomb phase) and a crystal with bro- 
ken lattice symmetries. Our high-precision Monte Carlo 
(MC) data allow independently to locate the transition 
starting from both phases. We determine the critical 
exponents and find that they are surprisingly close to 
those of a tricritical point. In a related recent work, 
Bergman et al. 01 further argue that the very existence 
of the Coulomb phase guarantees that this transition lies 
outside the standard LGW framework (see also Ref. 0). 
The Coulomb phase in turn owes its presence to the ab- 
sence of unpaired sites (monomers) 0- In this spirit, it is 
the suppression of such topological defects which opens 
the way to new types of critical behavior. Indeed, such 
"topological engineering" has p reviously been used in the 
context of liquid crystals [10|, and more recently for a 
non-linear a model [ll|; as we discuss before concluding, 
the numerical values of the critical exponents obtained 
in Ref. ^3 are not consistent with the ones reported here. 

The model is a 3d extension of the one studied in Ref.0 
on the square lattice. Configurations are dimer coverings 
of the simple cubic lattice of volume N = sites, with 
L the linear dimension. Dimers are hard-core and close- 
packed, i.e. every lattice site is part of one and only one 
dimer. Interactions favor aligment of nearest neighbours 
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(n.n.) dimcrs on plaquettes of the lattice: 

E=- ^ nj, +n//, (1) 

plaquettes 

with n|[,n= and n// denoting the number of plaquettes 
with parallel n.n. dimers in the x^y and z directions. 
Simulations (up to = 96"^) are performed with a recent 
MC directed- loop algorithm [l^ . 

At T = 0, the dimcrs order in columns to minimize 
the energy, resulting in a 6-fold degenerate ground state. 
The associated order parameter is a 3-component vector 
TO"(r) = {—y°"na{v)^ withnQ(r) — 1 for a dimcr pointing 
in direction a Cz x,y, z at site r, and otherwise. Naively, 
one would expect a high-T phase with (m) = and expo- 
nentially decaying dimer correlations. However, as shown 
by Huse et al. , at T = oo the system is in a "Coulomb 
phase" , with no true long-range order but with dipolar 
dimer-dimer correlations. To see this, the appropriate 
variable is the "electric field" Q Ea{r) = ( — )''?7Q(r). 
This field satisfies V • E = (— )'' = ±1, as the dimers 
are close-packed. The Coulomb phase can be character- 
ized in the continuum by an effective "electrostatic" 
action 8=^1 drE^ (r) which generates the dipolar cor- 
relations. Dimer fluxes (p = /j^E • dS through the planes 
perpendicular to the units vectors are conserved quanti- 
ties and vanish on average. One easily shows that flux 
fluctuations allow to calculate K 

{<t>')lL = ^m + {<t>l) + m) = l/K. (2) 

A close similarity with a 3d XY model can be seen 
through a duality transformation |0 in which V • E = 
(— is enforced by an angular Lagrange multiplier 6 at 
each site. The discrete sums on E^ are then performed 
by a Poisson formula, resulting in an XY interaction (Vil- 
lain form) between the variables. In this language, the 
Coulomb phase corresponds to an ordered phase with 
broken 0(2) symmetry for 6*, and \/K \s the associated 
spin stiffness. 

We first present thermodynamic results. Fig. (left 
panel) shows the behavior of the specific heat per site 
Cv/N as a function of T. Two close-by peaks are ob- 
served around T 1.52 and T ~ 1.67. The first peak 
is much broader and does not diverge with system size: 
since it is already present and almost converged on small 
lattices L < 16 (not shown), it cannot be associated to 
any long distance or critical behavior. The second peak 
is much more characteristic of a phase transition: it di- 
verges with L, with a power-law- like envelope typical of 
second order phase transitions (see top right panel) . Note 
that this peak is extremely hard to detect since it is ab- 
sent on small systems (L < 16) and also very sharp. We 
interpret this peak as the signature of the direct tran- 
sition (see below) between the Coulomb and columnar 
phase. Our best estimate for the temperature of its di- 
vergence is T^^ = 1.676(1). To determine the nature of 
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FIG. 1: Left panel: Specific heat per site Cv/N as function of 
temperature T for different system sizes L. Top right panel: 
Zoom on the second peak. Bottom right panel: Scaling of 
the specific heat at T^^ — 1.676 as a function of system size 
L. The dotted line denotes a power-law fit for the 4 largest 
systems L = 48, 64, 80, 96 (see text). 

the transition, we also considered energy histograms and 
the energy cumulant Q defined as 1 — {E'^) /3(£'^)^. 
No sign of double peak is detected in histograms and the 
energy cumulant saturates to 2/3 at the transition point: 
this indicates that the transition is not first order. 

Let us now consider the high-T phase. The left inset 
of Fig. 121 displays typical data for for a sample L = 
32. K is finite in the whole high-T range (with a value 
K{T = oo) = 5.12(1) in agreement with Ref. 0) and 
diverges below T ~ 1.6 (rough estimate from the plot, 
and inset of Fig.OJ. We expect for a second order phase 
transition a scaling form of the type 

i^-i = (02)/L = L-^/(Li/^(T-rj) (3) 

where z is a scale exponent, ly the correlation length ex- 
ponent, / a scaling function and Tc the critical temper- 
ature. Dimensional analysis of the Coulombian action 
gives z = 1 and therefore, at the curves of LK^^ in- 
tersect for all L and the derivative L '^^j, scales as L^/'^ . 
Numerically an accurate crossing point (see main panel 
of Fig. I2J) is indeed obtained for z = 1, which is also a 
good check of the second-order nature of the transition. 
An estimate = 1.6745(5) can be obtained from the 
crossing of the curves for the largest L. In the Coulomb 
phase, dimer-dimer correlations are expected to be dipo- 
lar 01, and this is found to be indeed true all along the 
high-T phase. The prefactor in the dipolar form of the 
correlation functions varies as l/i^T, and we have checked 
with high precision that the value of K obtained from 
flux fluctuations perfectly coincides with that from corre- 
lations. Monomer-monomer correlators are also available 
in the simulations la, Il2l and we find that test monomers 
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FIG. 2: Stiffness A'^^ multiplied by L (obtained from {(j)^), 
see Eq. \^ versus T near the transition, for different system 
sizes L. Left inset: Stiffness for the whole high-T range 
for L = 32. Right inset: Scaling of the derivative LdK~^ / dT 
versus L in log-log scale for the estimated critical temperature 
= 1.6745. The dotted line denotes a power-law fit for the 
4 largest systems. 
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FIG. 3: Crossing of the columnar Binder cumulant B as a 
function of T near the transition point for different system 
sizes L. The sample L = 16 is again out of scaling. Left inset: 
Columnar order parameter m and stiffness (multiplied 
by 8) in the whole T range for L = 32. Right inset: Cumulant 
derivative \dB/dT\ at T^"' = 1.67525 versus L in log-log scale. 
The dotted line is a power-law fit for the 4 largest L. 



are dcconfincd from T = oo down to Tc, confirming the 
Coulombian nature of the phase. 

To probe the nature of the low-T phase, we calculate 
the columnar order parameter 

w=— ||^TO(r)|| (4) 

r 

and its Binder cumulant B = 1 — (m^) /3 Irn^Y' . The 
left inset of Fig. |21 shows the expectation value (m) for 
a sample L = 32 (for illustration, is again shown). 
Columnar order is observed to set in at low T. Binder 
cumulants in the main panel admit a crossing point for 
systems with different L, leading to an estimate T^°' = 
1.67525(50). Assuming the standard scaling form B = 
fiL^/^XT - Tc)), the derivative dB/dT should scale as 
L^/" at Tc. 

The previous findings and the agreement between the 
various estimates of clearly indicate that the model 
displays a single second- order phase transition between 
Coulomb and columnar phases. A straightforward choice 
for a LGW theory would be to use m as an order pa- 
rameter. However, this fails as the Coulomb phase is not 
just a simple liquid where all correlations decay expo- 
nentially. Rather, it retains algebraic dimer correlations 
whose dipolar nature, crucially, does not lead to a peak 
in the structure factor anywhere in Fourier space (un- 
like the analogous situation in 2d 0, 0|). Indeed, the 
natural variable with algebraic correlations is a coarse- 
grained "electric" field E] however, this variable exhibits 
no long-range order in either phase. Instead, it is the fluc- 
tuations of E that distinguish the two phases. It is also 
instructive to examine the transition from the Coulomb 



side with the dual angles 9. It allows to map the dimer 
problem onto a model of interacting vortex loops with a 
long-range 1/r potential. While these loops are dilute in 
the Coulomb phase, they have to "proliferate" to repro- 
duce a low-temperature crystal phase with frozen dimer 
positions (thus highly fluctuating dual variables 6). In- 
tuitively, the restoration of the 0(2) symmetry would be 
through an inverted 3d XY transition, which is however 
incompatible with the critical exponents found numeri- 
cally (see below). The crucial difference with a simple 
0(2) spin model can be traced back to the background 
electric charges ±1 which couple to the vortex loops and 
presumably affect their proliferation. We note the simi- 
lar analysis of Ref. 0, where an unconventional non-LGW 
transition is predicted in a closely related model. 

We now come to the universality class of the transition. 
The correlation length exponent v can be extracted from 
the scaling with L of stiffness LdK~^ /dT or Binder cu- 
mulant dB/dT derivatives at the critical temperature T^, 
which can be calculated thermodynamically in the MC 
process. Taking into account only the largest L > 48 (see 
insets of Figs. [21and|3J), we obtain compatible estimates 
= 0.50(4) and v'^°'^ = 0.51(3) (error bars take into ac- 
count stability of fits toward inclusion of smaller samples 
and uncertainty on Tc). The specific heat critical expo- 
nent a can be extracted from its scaling at the critical 
point: Cv(Tc)/N = cq + AL"'^'^ where A is a constant and 
Co the regular part at the transition (co is non negligible 
as can be seen for the L = 16 sample in Fig. A fit for 
the largest L (see inset in Fig. ^ gives a/v = 1.11(15), 
leading to a = 0.56(7). Hyperscaling a = 2 ~ i/d is thus 
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satisfied within error bars. The last independent expo- 
nent can be obtained from the scahng of the columnar 
order parameter at criticality m{Tc) ^ L~^/'^ or from the 
associated susceptibility x = (m'^) — {m,)'^ ^ L^/^ . Us- 
ing standard relations between critical exponents, wc ob- 
tain for the correlation function exponent rj = —0.02(5). 
This set of exponents excludes some simple 3d univer- 
sality classes (such as 0(2), 0(3) with or without cubic 
anisotropy) but are compatible with the universality class 
of an 0(n) tricritical point at its upper critical dimension 
d = 3 for which v = a = 1/2 and ij = (up to logarith- 
mic corrections). We also note that the value at Tc of 
the cumulant of the 3-dimensional order parameter m is 
compatible within error bars (see Fig. ^ with the value 
0.56982 ... of a tricritical 0(3) theory (for d > 3) [H. 
With dimers, we do not have direct access to the XY 
order parameter n(r) for the dual angles 9. We can how- 
ever investigate the fluctuations of the electric flux. 0^ is 
the integer-valued Noether charge associated to the 0{2) 
symmetry ( "total angular momentum" if z is interpreted 
as the time direction): cpz = n J dxdy{nidzn2 — n2dzni) 
where k appears in the "kinetic" term ■|(c?^n)^ of the 
dual 0(2) action. Above Tc the typical flux scales as 
Vl and the ratios {(j^t) / {4>l)'^ and {(f'^) / {(j)^)'^ are re- 
spectively equal to 3 and 5/3, in agreement with the 
Gaussian and 0(3)-symmetric nature of fluctuations in 
the Coulomb phase. We believe that the distribution 
of (f) is universal at Tc and we measured (0^) = 0.28(2) 
and {(f)^) = 0.25(4). The smallness of these quantities 
at the critical point means that the discrete nature of 
4> = {4>x, 4>yj4'z) cannot be neglected there. 

While it may be seen as the only way to reconcile the 
numerical results with a LGW analysis, the tricritical 
universality class is rather unexpected here as it would 
imply a "hidden" fine-tuning of parameters of the effec- 
tive action. It is also quite possible that the exponents 
found are close, but not equal, to tricritical exponents, 
thereby defining a new "non-LGW" 0, 0| universality 
class. The absence of monomers is in fact very simi- 
lar to the absence of "hedgeho gs" in the models studied 



by Motrunich and Vishwanath which also display a 
transition from a broken-symmetry phase to a Coulomb 
liquid. However, our critical exponents do not match 
those of Ref.0. This discrepancy might be due to one or 
several factors: i) lattice cubic anisotropies (not present 
in Ref. ^3 but potentially relevant) exist in our model 
which admits 6 ground-states related by cubic symmetry, 
ii) the simulations of Ref. may not be in the scaling 
regime, iii) the proximity of a tricritical point (for in- 
stance at finite monomer doping), could affect finite-size 
estimations of the exponents and hide the true critical be- 
havior. Finally, even though we simulate systems with N 
up to 96^, a first-order transition with a large correlation 
length can never be excluded from finite-size simulations. 

In conclusion, our study has established the crystal- 
lization in the cubic dimcr model as an example of 



an unconventional phase transition in a classical model 
in 3 dimensions. The conventional LGW approach to 
phase transitions is currently under attack from many 
sides 0, 0, 0, 0| , but shows considerable resilience in 
microscopic models 0. From this perspective, our re- 
sults are promising since they cannot be described by a 
strict application of the LGW scheme of an expansion 
in terms of the low-T order parameter. Further analyti- 
cal calculations and numerical tests are needed to extend 
this study, notably to investigate the possible proxim- 
ity to a tricritical universality class. This will hopefully 
allow a deeper understanding of the limits of a central 
concept in statistical physics. 
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